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Abstract 

We study the Whitham equations for the Camassa-Holm equation. The equations 
are neither strictly hyperbolic nor genuinely nonlinear. We are interested in the 
initial value problem of the Whitham equations. When the initial values are given 
by a step function, the Whitham solution is self-similar. When the initial values are 
given by a smooth function, the Whitham solution exists within a cusp in the x-t 
plane. On the boundary of the cusp, the Whitham solution matches the Burgers 
solution, which exists outside the cusp. 
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1 Introduction 

The Camassa-Holm equation 

Mt + (3m + 2z/)m^ - e^(n3;3;t + 2m^m^.^. MMxxa;) = , M(x,0;e) = Uo(a;) (1) 

describes waves in shallow water when surface tension is present [2]. Here, 
is a constant parameter. The solution of the initial value problem ([T| will 
develop singularities in a finite time if and only if some portion of the positive 
part of the initial "momentum" density uo{x) — e^UQ^x) + v lies to the left of 
some portion of its negative part [10]. In particular, a unique global solution 
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is guarenteed if uo{x) — e'^u^^x) + v does not change its sign. These are the 
non-breaking initial data that we are interested in throughout this paper. 

Ahhough the zero dispersion hmit of equation ([T]) has not been estabhshed, 
some of its modulation equations (i.e., Whitham equations) have been derived. 
The zero phase Whitham equation is 

Ui + (3m + 2v)ux = , u{x, 0) = uo{x) , (2) 

which can be obtained from ([T| by formally setting e = 0. 

The single phase Whitham equations have been found in [T] and they can be 
written in the Riemann invariant form 

Uix + Xi{ui,U2,U3)Uix = ioT - U < U3 < U2 < Ui , (3) 

where ^ 

\i{Ui,U2,U3) = Ui+U2+Us + 2u - -— , (4) 

and 

■t] + u 



I{ui,U2,U3)= dr] . 

y{v + '^)iui-r]){u2-r]){r]-U3) 

The constraint — z/ < M3 < -^2 < is consistent with the non-breaking initial 
data mentioned in the first paragraph of this section. The integral I can be 
rewritten as a contour integral. Hence, it satisfies the Euler-Poisson-Darboux 
equations 

d'^I dl dl 

2{ui-Uj) =- — , z,j = 1,2,3 (5) 

OUiOUj OUi OUj 

since the integrand does so for each 7] ^ Ui. The Hamiltonian structure of 
the single phase Whitham equations for the Camassa-Holm equation was also 
obtained in [T] in terms of Abelian integrals. The higher phase Whitham equa- 
tions can also be derived using this structure. 

In this paper we will study the evolution of the Whitham solution from the 
zero phase to the single phase. 

This problem is similar to that of the zero dispersion limit of the KdV equation 
[T1I81IT71 

ut + 6uux + e^Uj^xx = , u{x,0;e) = uo{x) . (6) 
It is known that the zero phase Whitham equation for the KdV equation is 

ut + 6uux = , (7) 

which is equivalent to ^ for the Camassa-Holm equation. The single phase 
Whitham equations for the KdV equation are [3|HIT8] 



Uix + IJ'i{UuU2,U3)Uix = , for U3<U2<Ui, 
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where 



/ 

fli{Ui, U2, U3) = 2[Ui+U2 + U3- 



and 



/(mi,M2,M3) = / I, _ _ ^ dr] 

Juz 



These equations are also similar to ^ and Q for the Camassa-Holm equation. 

In the KdV case, the evolution from the zero phase to the single phase has 
been studied in [H]. There, the Euler-Poisson-Darboux equations ([s]) have 
played an important role. The same equations have also played a crucial role 
in the study of the transition from the single phase to the double phase in [3] . 

Although both the Camassa-Holm equation ([T| and the KdV equation ^ 
are dispersive approximations to the Burgers equation (|2]) or ([T]), there are 
significant differences in the limiting dynamics. The biggest difference is that 
the Whitham equations ([s]) for the former equation are non-strictly hyperbolic 



(cf. (32)) while the Whitham equations (|8]) for the latter equation are strictly 
hyperbolic [9]. Non-Strictly hyperbolic Whitham equations have also been 
found in the higher order KdV flows PT|12II13] and the higher order defocusing 
NLS flows Self-similar solutions of these Whitham equations have been 
constructed. They are remarkably different from the self-similar solutions of 
the KdV- Whitham equations pT]|T2] or the NLS- Whitham equations [6], both 
of which are strictly hyperbolic. 



In this paper, we will modify the method of paper p!l] so that it can be used 
to solve the non-strictly hyperbolic Whitham equations ^ when the initial 
function is a smooth function. We will then study the evolution from the zero 
phase to the single phase for smooth initial data. When the initial function 
is a step-like function, we will use the method of paper Plll|12j to study the 
same evolution. 



The organization of the paper is as follows. In Section 2, we will introduce 
an initial value problem which describes the evolution of phases. We will also 
discuss how to use the hodograph transform to solve non-strictly hyperbolic 
Whitham equations. In Section 3, we will study the properties of the eigen- 
speeds of the single phase Whitham equations. We will study the initial value 
problem when the initial function Uo{x) is a step-like function in Section 4 and 
when Uq{x) is a smooth decreasing function in Section 5. 
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2 An Initial Value Problem 



We describe the initial value problem for the Burgers equation ^ and Whitham 
equations ^ as follows (see Figure 1.). Consider a horizontal motion of the 
initial curve u = uq{x). At the beginning, the curve evolves according to the 
Burgers equation (|2]). The Burgers solution breaks down in a finite time. Im- 
mediately after the breaking, the curve develops three branches. Denote these 
three branches by Ui, U2, and u^. Their motion is governed by the Whitham 
equations ([s]). As time goes on, the Whitham solution may develop singular- 
ities and more branches are created. However, our focus is on the evolution 
of the solution of the Whitham equations from the one branch regime to the 
three branch regime. 





u 










u 








X.(t) X^t) X 



Fig. 1. Profile of the Burgers and Whitham solutions. The Burgers solution u of (2) 
lives in the single value regions while the Whitham solution ui, U2 and of (3) 
reside in the multiple value region. 

The Burgers solution m of ([2]) and the Whitham solution -ui, U2, U3 of ^ must 
match on the trailing edge x = X-{t) and leading edge x = x^{t). We see from 
Figure 1 that 



Ui = u 
U2 = M3 



(9) 



must be imposed on the traihng edge, and that 



Ui = U2 
U3 = U 

must be satisfied on the leading edge. 



(10) 
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In this paper, we consider the initial function uq{x) that is monotone. Since the 
Burgers solution will not develop any shock if Uq{x) is an increasing function, 
we will focus on decreasing initial functions. Denoting the inverse function of 
uo{x) by f{u), the Burger equation ^ can be solved using the method of 
characteristics; its solution is given implicitly by a hodograph transform, 

x= {3u + 2u)t + f{u) . (11) 



The solution method (11) has been generalized to solve the first order quasi- 
linear hyperbolic equations which can be written in Riemann invariant form 
and which are strictly hyperbolic 



duj 

In 



+ sAui, 



dui 



1,2, ■■ ■ ,n. 



(12) 



The strict hyperbolicity means that the wave propagation speeds Sj's do not 
coincide. 

We include Tsarev's theorem for completeness p^lllG] . 
Theorem 1 If Wi{ui,U2, ■ ■ ■ ,Un) 's solve the linear equations 



dwj 



with 



Aij{wi 



dSj_ 

dui 



for i, j = 1,2, ■ ■ ■ ,n and i ^ j , then the solution {ui{x, t), 
hodograph transform, 



(13) 
(14) 

, Un{x, t)) of the 



X = Si{Ui, U2,--- , Un)t + Wi{Ui, U2, ■ ■ ■ , Un) (15) 

satisfies equations p^. Conversely, any solution {ui{x,t), ■ ■ ■ ,Un{x,t)) of 



equations (12) can he obtained in this way in the neighborhood of (a;o,to) 
which are not all vanishing. 



The strict hyperbolicity, i.e., Sj 7^ Sj for i 7^ j, of (12) is assumed to ensure 
that Ay's of (14) are not singular. 

The result is classical when n = 2. 

The validity of this theorem hinges on two factors. First, the linear equations 
(13) must have solutions. Secondly, the hodograph transform (15) must not 
be degenerate, i.e., it can be solved for Ui S cLS functions of x and t. One 
interesting observation is that the Jacobian matrix of ( 15 ) is always diagonal 
on the solution {ui{x, t), - ■ ■ ,Un{x,t)). 
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Corollary 2 At the solution {ui{x,t), ■ ■ ■ ,Un{x,t)) of x = Si{ui,U2, ■,Un)t + 
Wi{ui,U2, ■ ■ ■ , Un), i = 1,2, ■ ■ ■ ,n, the partial derivatives 

d{sit + Wj) _ ^ 
duj 

fori 7^ j. 
Proof. 



d{sit + Wj) _ dsi^^ ^ dwj 
duj duj duj 



Aj [{Sit + Wi) - {Sjt + Wj)] 

. 



□ 



Another aspect of Theorem [T] is that it is a local result. Solutions produced 
by the hodograph transform are, in general, local in nature. However, global 
solutions can still be obtained if the conditions of the theorem are satisfied 
globally pii] . 

The Whitham equations ^ will be shown to be non-strictly hyperbolic, i.e., 
Aj's coincide at some points {ui,U2, M3) where 0<U3 + h'<U2 + i^<ui + h'. 
However, Theorem [l] can still be applied to equations ^ since the functions 
Aij^s of (14) are still non-singular for the Whitham equations (|3|, even at the 



points of non-strict hyperbolicity. 
Lemma 3 

' Aj - Aj 2 {Xj - 'y){ui - Uj) 
where 'j = ui + U2 + U3 + 2z/. 
Proof. 

By Q, we calculate 



dXj ^ 2{u,-u,)dl,J ^ Idl^J 
duj duj {dujy 

where we have used ([s]). Hence, we get 
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_ 1 {Xj - 7) - - Uj) 
where we have used Q to express I /duj- n 



3 The Single Phase Whitham Equations 



In this section, we will summarize some of the properties of the speeds Aj's of 
Q for later use. 



Function / of Q is a complete elliptic integral; indeed, 

r( ^ 2(-»3 + Z/)n(p, S) 

where 11 (p, s) is the complete integral of third kind, and 

U2 + V ' (Mi - Uz){u2 + V) 



P 



(16) 



(17) 



Properties of complete elliptic integrals of the first, second and third kind are 
listed in Appendix A. 



Using the well known derivative formulae (A. 3) and (A. 4), one is able to 
rewrite Aj of ([i]) as [1] 



Ai(mi, M2, U3) = ui + U2 + u-i + 2v + 2{ui - U2) 



(-»3 + y)Ti{p, s) 
{u2 + v)E{s) 



A2(Mi,M2,%)=Mi + U2 + U3 + 2z/ + 2(m3-M2)7^7^^ f^^^'^wL ^ , (18) 



A3(mi, M2, Ms) = Ml + M2 + Ms + 2z/ + 2(-U2 - Ms) 



E{s) - {l-s)K{s) 

(m3 + z/)n(p, s) 



{u2 + u)[E{s)-K{s)] ■ 
Here K{s) and E{s) are complete elliptic integrals of the first and second kind. 



Using inequalities (A. 9), we obtain 
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Ai - (mi + M2 + Ms + 2z/) > 
A3 - {Ui + U2 + M3 + 2z/) < 



(19) 
(20) 
(21) 



for ui > U2 > U3 > —V. In view of (A. 5 A.8) and (A. 10 A. 11), we find that 
Ai, A2 and A3 have behavior 



(1) At U2 = U3, 

Ai(ui, U2, Us) = 3ui + 2z/ , 

A2(ui, U2, Us) = Xsim, U2, Us) = Ml + 2m3 + 2u - 



(22) 



(2) At Ml = M2, 



Lemma 4 



Ai(Mi, M2, M3) = A2(Mi, M2, M3) = 2Mi + M3 + 2i/ 
A3(Mi,M2,M3) = 3m3 + 2u . 



dXs ^ 3 A2 - A3 ^ 0X2 

dus 2 M2 - M3 du2 



for < M3 + z/ < M2 + < Ml + z/. 
Proof. 

Comparing formulae ^ and (18), we use (16) to obtain 

dl \Jiui - M3)(m2 + U) 



dU2 {U2 - Us){Ui - U2) 

Differentiating ^ for A2 and using (24) yields 



\E-{1- s)K] 



(23) 



(24) 



dX, 



du2 (duj)^ 

{us + u)U 
iu2 + iy)[E-{l~s)K]^ 

X ([. + - 2'^^][E - (1 - s)K] + s{l - s)k] . 

I Ml - M3 Ml - M3 J 



(25) 



Using formulae (18) for A2 and A3, we obtain 

2(m2 - M3)(m3 + Z/)nM(Mi, M2, M3) 



A2(Mi, M2, M3) - A3(Mi, M2, M3) 



{u2 + u){K -E)[E-{l-s)K] 



(26) 



8 



where 



M(Mi,M2,M3) = [IH \E - [ 



[1 - s)]K . (27) 



We then obtain from (25) and (26) that 



0X2 3 A2 — A3 



(M3 + Z/)(4-3s)nE2 



du2 2 U2 — U3 



{u2 + u)[K -E][E-{l-s)K]n^^ ""^^E^^ 



E' 4 -3s 



> 



where the inequahty follows from the negativity of the function in the bracket 
(c.f. (4.18) of [S]). This proves part of the lemma. The other part can be 
shown in the same way. 



□ 



We conclude this section with a few calculations. We use (A.l) and (A. 2) to 



calculate the derivative of (27) 



gM(Mi,M2,M3) 
dU2 



d^M{ui,U2,U3) 

dUn 



2(mi - U3) {u2 + 

U2 + l^ 

(ui + u)[A{ui + z/) - 2(m2 + + jus + u)]E 

4(Mi - %)(«! - U2){U2 + VY 



(28) 
(29) 



Finally, we use the expansions (A. 5 A. 6) for K and E to obtain 



N 7rJ(M2-M3)s 1 3(ni-n2). 2 
= 2 1 2(.7^ + 16 



Ui - U3 



Ui - U3 



TT 

2 



+-(1 

128^ 

[ ^^ + '^ +-(1- 
^2{ui + u) 16 ^ 

+4(1 



3{ui - U2] 

Ui - U3 



128 



Ui - U3 



(30) 



where we have used formula (17) for s in the last equality. 
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4 Step-like Initial Data 



In this section, we will consider the step-like initial data 



a X < 
h X > 



(31) 



for equation ([T]). Since the solution of ([2]) will never develop a shock when 
a < 6, we will be interested only in the case a > b. We classify the initial data 



(31 ) into two types: 



• (I) a + z/ > 4(6 + z/) , 

• (II) a + z/<4(6 + z/). 

We will solve the initial value problem for the Whitham equations for these 
two types of initial data. 



4.1 Type I: a + V > A{h + v) 



2.0 



1 .5 



u 



1.0 



0.5 



0.0 













u 
























x/t=a 








\u 


3 


u 













x/t 

Fig. 2. Self-Similar solution of the Whitham equations for a = 2, 6 = and v = 1 /20 
of type I. 



Theorem 5 (see Figure 2.) For the step-like initial data uq{x) of (31) with 



< 6 + z/ < (a + z/)/4, the solution of the Whitham equations is given by 

ui = a , a; = A2 (0,^2,^3) t , a; = A3(a, M2, U3) ^ (32) 
for (3a — z^)/4 < x/t < a and by 

Ui = a , X = A2(a, U2, b) t , U3 = b (33) 
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for at < X < (2a + b + 2z/)t, where a = X2{a,u*,b) and u* is the unique 
solution U2 of X2{a,U2,b) = \^{a,U2,b) in the interval b <U2 < a. Outside the 
region (3a — i/)/4 < x/t < 2a + 6 + 2u, the solution of the Burgers equation 
is given by 

u = a x/t < {3a- u) /A (34) 

and 

u = b x/t>2a + b + 2u . (35) 

The boundaries x/t = (3a — i/)/4 and x/t = 2a + b + 2u are the traihng 
and leading edges, respectively, of the dispersive shock. They separate the 
solution into the region governed by the single phase Whitham equations and 
the region governed by the Burgers equation. 

The proof of Theorem [5] is based on a series of lemmas. 



We first show that the solutions defined by formulae (32) and (33) indeed 
satisfy the Whitham equations ^ [16] . 



Lemma 6 (i) The functions Ui, U2 and determined by equations (32) 

give a solution of the Whitham equations ^ as long as U2 and us can be 

solved from (32) as functions of x andt. 
(a) The functions Ui, U2 and determined by equations (33) give a solution 

of the Whitham equations as long as U2 can be solved from (33) as a 

function of x andt. 

Proof. 

(i) Ml obviously satisfies the first equation of (|3|. To verify the second and 
third equations, we observe that 

^ = ^ = (36) 
du-i du2 

on the solution of (32). This follows from Lemma [3] 



We then calculate the partial derivatives of the second equation of (32) with 
respect to x and t. 

9X2 / N 9X2 , ^ 

dU2 OU2 

which give the second equation of (|3|. 
The third equation of (|3| can be verified in the same way. 
(ii) The second part of Lemma [6] can easily be proved in a similar manner. 
□ 
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We now determine the trailing edge. Eliminating x and t from the last two 

Is 

A2(a, ^2,^3) - A3(a,U2,M3) = . (37) 



equations of (32) yields 



In view of formula (26), we replace (37) by 



M(a, M2,M3) 



. 



(3J 



Therefore, at the trailing edge where U2 = M3, i.e., s = 0, equation (38), in 

3(a - M3) 



view of the expansion ( 30 ) , becomes 

1 



^^ + ^ +-[1 
2(a + z/) 16 ^ 



which gives U2 = = {a — 3z/)/4. 



a — M3 







Lemma 7 Equation (38) has a unique solution satisfying U2 = U3. The solu- 



tion is M2 = M3 = (a — 3iy)/A. The rest of equations at the trailing edge 
are Ui = a and x/t = A2(a, (a — 3z/)/4, (a — 3z/)/4) = (3a — z^)/4. 



Having located the trailing edge, we now solve equations (32) in the neigh- 



borhood of the trailing edge. We first consider equation (38). We use (30) to 



differentiate M/s at the trailing edge ui = a, U2 = U3 = {a — 3z/)/4, to find 



a M d M 
du2 du^ 



TT 

ia + v) 



which shows that equation (|38|) or equivalently ([37j) can be inverted to give 

(39) 



M3 as a decreasing function of U2 

U3 = A{U2) 

in a neighborhood of U2 = = {a — 3z/)/4. 



We now extend the solution (39) of equation (37) in the region a > U2 > 



(a — 3i/)/4 > M3 > 6 as far as possible. We deduce from Lemma ^ that 



d\2 



> 



dU3 



< 



(40) 



on the solution of (37). Because of (36) and (40), solution (39) of equation 



(37) can be extended as long as a > M2 > (a — 3z/)/4 > M3 > 0. 



There are two possibilities: (i) U2 touches a before or simultaneously as U3 
reaches h and (ii) M3 touches h before U2 reaches a. 
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It follows from (23) that 

A2(a, a, Us) > Xs^a, a, for 6 < M3 < a . 

This shows that (i) is impossible. Hence, M3 will touch b before U2 reaches a. 
When this happens, equation (37) becomes 

X2{a,U2,b) - X3{a,U2,b) = . (41) 



Lemma 8 Equation (4I) has a simple zero in the region b < U2 < a, counting 
multiplicities. Denoting the zero by u* , then A2(a, U2, 6) — A3(a, U2, b) is positive 
for U2 > u* and negative for U2 < u* . 



Proof. We use (26) and (29) to prove the lemma. In equation (26), K — E 
and E — [1 — s)K are all positive for < s < 1 in view of (A. 9). We claim 
that 

M(a, U2, 6) = for M2 = 6 , M(a, U2, 6) < for U2 near b , 



and 



M(a, ^2, &) > for ^2 = a • 



The equality and the first inequality follow from expansion (30) and a + z/ > 
4(6 + v). The second inequality is obtained by applying (A.7) and (A. 8) to 
(0. 

We conclude from the two inequalities that M(a, U2, b) has a zero uib < U2 < a. 
This zero is unique because M{a,U2,b), in view of (29), is a convex function 
of U2. This zero is exactly u* and the rest of the theorem is proven easily. □ 



Having solved equation (37) for as a decreasing function of U2 for (a — 
3i^)/4 < M2 < u*, we turn to equations (32). Because of (36) and (40), the 
second equation of (32) gives U2 as an increasing function of x/t, for (3a — 
z/)/4 < x/t < a, where 

a = A2(a, u*, b). 

Consequently, is a decreasing function of x/t in the same interval. 



Lemma 9 The last two equations of (32) can be inverted to give U2 and M3 as 
increasing and decreasing functions, respectively, of the self- similarity variable 
x/t in the interval (3a — u) < x/t < a, where a = X2{a,u*,b) and u* is given 
in Lemma\^ 



We now turn to equations (33). We want to solve the second equation when 



x/t > a or equivalently when U2 > u*. According to Lemma [8| X2{a,U2,b) 
X3{a,U2,b) > for u* < U2 < a, which, by Lemma |4| shows that 



dX2{a,U2,b) 
duo 



> . 
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Hence, the second equation of (33) can be solved for U2 as an increasing func- 
tion of x/t as long as m* < M2 < a. When U2 reaches a, we have 



x/t = \2{a, cL,b) = 2a + 6 + 2z/ 



where we have used ([23j) in the last equality. We have therefore proved the 
following result. 



Lemma 10 The second equation of (33) can be inverted to give U2 as an 



increasing function of x/t in the interval a<x/t<2a + b + 2iy. 
We are ready to conclude the proof of Theorem [5] 



The Burgers solutions (34) and (35) are trivial. 



According to Lemma [oj the last two equations of (32) determine U2 and us as 
functions of a;/t in the region (3a — f )/4 < x/t < a. By the first part of Lemma 
6| the resulting ui, U2 and U3 satisfy the Whitham equations ([3]). Furthermore, 
the boundary condition ^ is satisfied at the trailing edge x/t = (3a — v)/A. 



Similarly, by Lemma [10| the second equation of (33) determines U2 as a func- 



tion of x/t in the region a < x/t < 2a + b + 2u. It then follows from the 



second part of Lemma 6 that Ui, U2 and % of (33) satisfy the Whitham equa- 



tions (|3|. They also satisfy the boundary condition (10) at the leading edge 
(2a + 6 + 2v)t. 



X 



We have therefore completed the proof of Theorem |5] 



A graph of the Whitham solution for the initial data (31 ) of type I is given in 



Figure 2. It is obtained by plotting the exact solutions of (32) and (33) 



12 Type II: a + V < A(b + u) 



Theorem 11 (see Figure 3.) For the step-like initial data (31) with < {a-\- 
z/)/4 < 6 + z/ < a + z/, the solution of the Whitham equations is given by 

ui = a , X = A2(a, U2, b) t , = b (42) 

for A2(a, b,b) < x/t < A2(a, a, b), where A2(a, b,b) = a-\-2b + 21/ — 4(a — b){b-\- 
v)/{a + v) and A2(a, a, 6) =2a + b + 2v. Outside this interval, the solution of 
is given by 

u = a x/t < X2{a,b,b) 

and 

u = b x/t > \2{a, a, b) . 
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Fig. 3. Self-Similar solution of the Whitham equations for a = 2, b = 0.8 and 
1/ = 1/20 of type II. 



Proof. 



We will give a brief proof, since the arguments are, more or less, similar to 
those in the proof of Theorem [5j 



It suffices to show that A2(a, U2, b) is an increasing function of U2 for b < U2 < a. 



Using the inequality (A. 9) to estimate the right hand side of (28), we obtain 



dM{a, U2, b) 
duo 



> 



{a + u){u2 - b)E{s) 
2{2 -s){a- by {u2 + uy 



{2{u2 + i^) + 2{b + u)- 



-{a + u)} > 



for b < U2 < a, where we have used (a + z/)/4 < 6 + 1/ in the second inequality. 
Since M(a, U2, 6) = at M2 = & in view of (30), this implies that M(a, U2, b) > 



for 6 < ^2 < a. It then follows from (26) that A2(a, U2, b) — \^{a, U2, b) > 0. By 
Lemma |4| we conclude that 



dX2{a,U2,b) 
duo 



> 



for b < U2 < a. 



□ 



A graph of the Whitham solution for initial data (31) of type II is given in 



Figure 3. It is obtained by plotting the exact solution of (42). 
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5 Smooth Initial Data 



In this section, we will study the initial value problem of the Whitham equa- 
tions when the initial values are given by a smooth monotone function Uq{x). 
Since the Burgers solution of ^ will never develop a shock when uo{x) is 
an increasing function, we will be interested only in the case that uo{x) is a 
decreasing function. 

We consider the initial function uo{x) which is a decreasing function and is 
bounded at a; = ±00 



lim Uo{x) 



lim Uq(x) = b 



(43) 



By Theorem [T] and Lemma [3} we can use the hodograph transform, 
X = Xi{ui, U2, uz)t + Wiiui, U2, U3) , i = 1, 2, 3, 



(44) 



to solve the Whitham equations (|3|). Here, Wj's satisfy a linear over-determined 



system of type ( 13 ) 



dwj 

dUn 



Bijiwi -Wj) , 



(45) 



where -Bj/s are given in Lemma |3 



The boundary conditions on U7j's are obtained by observing that the hodograph 
solution (44) of the Whitham equations (|3| must match the characteristic 
solution (11) of the Burgers equation (|2|) at the trailing and leading edges in 



the fashion of ([9 10). By (22 23), Wi's must satisfy the boundary conditions 



Wi{Ui,Ui,Us) 
Wi{Ui,U3,U3) 



W2{Ui,Ui,U3), 
/(«3), 
^ /(«l), 

-- 1^3(^1, M3, Us), 



(46) 
(47) 



where f{u) is the inverse of the initial function Uq{x). 

Analogous to the KdV case p3p5] . equations ([45]) subject to boundary con- 



ditions (46 47) are related to a boundary value problem of a linear over- 



determined system of Euler-Poisson-Darboux type (cf. (p 
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2{ui - Uj) 



d'^q dq dq 



duiduj dui duj 
q{u,u,u) = f{u) , 



hj = 1,2,3, 



(48) 
(49) 



for i,j = 1, 2, 3. The solution is unique and symmetric in ui, U2,U3. It is given 
explicitly by [1] 



q{ui,U2,U3) 



Theorem 12 // g(Mi, ^2, Ms) 2S a solution of (48) and (49), then 
Wi{ui,U2,u^)'s given by 



Wi = [Xi{Ui,U2,U3) - 'j\ — hg(Ml,M2,M3j 



(51) 



where 'j = Ui + U2 + u-^ + 2v, solve equations (45) and satisfy boundary con 



ditions (46 4"^) 
Proof. 



By Lemma and (48), we obtain 



[A, - 7] 



dq 

dui 



dq 



' dq 



dq 

dui 



+ [A. - 7] 



Using (51), we calculate 



d'^q 

duiduj 



(52) 



dwi dXi dq d'^q dq dq 

duj duj dui * 



duiduj duj dui 



dq dq dq 



dui 



Substituting these into (52), we find that Wj's satisfy (45). 



Finally, we shall check the boundary conditions (46 47). We only consider the 



leading edge, and the trailing edge can be handled in the same way. 



Since q{ui,U2,u-i) is symmetric in ui, U2 and the first condition of (46) 



follows from (51) and (23). 
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For the second condition, it follows from (23) and (51) again that 



, dq 

W3{Ui,Ui,U3) = 2{U3 - Ml)- h q . 

0U3 

Differentiating this with respect to Ui yields, 



(53) 



dui 



dq , d'^q 

-2^ + 2(M3-ni)[ 



dm 



+ 



d'^q 



dq dq 



duidu^ du2dus dui du 



, 



where we have used (48) in the last equality. Since w^iui^Ui, u^) is independent 



of (46). 



of Ml, we replace Mi by M3 in (53) and use (49) to obtain the second condition 



□ 



Theorem 12 has been reported in [I] in the case of z/ = 0. 



In the rest of this section, we study the hodograph transform (44) with Wj's 



given by (50) and (51). We shall show that the transform can be solved for 
Ml, M2 and M3 as functions of (x,t) within a cusp in the x-t plane. 

Since tj, = — [3 min (mo(x))]~^ is the breaking time of the Burgers solution of 
(|2|, the breaking is caused by an inflection point in the initial data. If Xq is 
this inflection point, then {xh,tb) is the breaking point on the evolving curve 
where Xb = xo+ [3mo(xo) + 2z/]tb, and U is the breaking time. Without loss of 
generahty, we may assume Xb = 0, tb = and denote mo(0) by u. The effect 
of these choices is that we are starting at the breaking time, and the evolving 
curve is about to turn over at the point (0, u) in the x-u plane. It immediately 
follows that 

/(m) = fiu) = f\u) = , (54) 

where x = f{u) is the inverse function of the decreasing initial data u = Uo{x). 
On the assumption that x = f{u) has only one inflection point, it follows from 
the monotonicity of the function /(m) that 



f'iu) 



< M > M 
= U = U 
> M < M 



(55) 



Under a little bit stronger condition than (55), we will be able to show that 



hodograph transform (44) can be inverted to give Mi, M2 and M3 as functions 



of (x, t) in some domain of the x-t plane. 



Theorem 13 Suppose Uq{x) is a decreasing function satisfying (43) with a 



u > b + u > 0. If, in addition to (54), the inverse function f{u) satisfies f' 



M 
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< for h < u < a, then transform (44) with Wi 's given by (50) and (51) can 
be solved for ui, U2 and as functions of (x, t) within a cusp in the x-t plane 
for all t > 0. Furthermore, these Ui, U2 and U3 satisfy boundary conditions 



(46 4"^) on the boundary of the cusp. 



The proof is based on a series of lemmas. The organization is as follows: we 



eliminate x from transform (44) to obtain two equations involving ui, U2, 



U3, and t. These two equations can be shown, for each fixed time after the 
breaking, to determine Ui and as decreasing functions of U2 within an 
interval whose end points depend on t. Substituting ui and as functions of 
U2 into the hodograph transform, we find that, within a cusp in the x-t plane, 
U2 is a function of {x,t), and so, therefore, are ui and U3. 



First, we conclude from formula (50) 



given by (50) satisfies 



Lemma 14 If f{u) satisfies the conditions of Theorem 13, then q{ui,U2,Us) 

<0, t,j,k = 1,2,3. 



duidujduk 



Eliminating x from (44) yields 



(Alt + - (A2t + W2)=0 

{\2t + W2) - {\zt + W^)={] 



(56) 
(57) 



Using (51) for wi and W2, and (18) for Ai and A2 , we write 



(Alt + w^) - (Aat + W2) = (Ai -l)\t+^]~ (A2 ~l)[t+ p-] 

OUi OU2 



where 



2{ui - U2){U3 + u)Il 
{U2 + v)E 

sE U2 + u 



F(ui,U2,U^) 

\t ' 



dui E — (1 — s)K ui + u du2 



Similarly, we use (51) for 'W2 and to write 



(5J 



da da do 

{\2t + W2) - (Aat + w^) = (A2 - A3)(t +Q^) + (^3 - - 



2{U2 - M3)(^3 + 

{u2 + u)[K-E][E-{l-s)K] 



G{Ui,U2,U^) 
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where 



G = Miu,,U2,Us)it + 1^) - 2(^2 - Us)[E - (1 - s)K]-^ . (59) 

OU2 0U20U^ 



In the derivation, we have used formula (26) for A2 — A3, formula (18) for A 



and equation (48) 



Since (A.9) implies that -ft'(s)—i?(s) > and £'(s) — (1 — s)-ft'(s) > 0, equations 

(60) 



(56) and (57) are equivalent to 

F(ui, M2, M3) = , ^(Mi, M2, %) = 



for < s < 1. 



5.1 The trailing edge 



We first study the trailing edge. We use (A. 5), (A. 6), and (30) to expand 



dui {ui + v) ^ du2 4(mi + z/) ^ du2 

3(m2 + ^) , \ „2 , r^i 3n 



and 



32(mi + v] 



OUo 



(61) 



G=- 



TT I ^ M2 + Z/ _^ ^ 3(M^^2^)j(^^ . 



2 [ 2{ui + u) 16 ■ Ml — M3 



dU2 



Ui + V 



du2du3 



— (l-5^^i^)(t+|^^ 
2\l28^ ui-m'^ du2' 

{ui - U3){u2 + v) d^q \ ^ 
8(ui + z/) du2du3 J 



0(s^ 



(62) 



Taking the limits of F = and G'/s^ = Oass— *0 and simplifying the results 
a bit, we obtain the equations governing the trailing edge 

U{u,,u,) := {u,+u) it+ )+2(n3+.) (t+ ^^("l; ) = , (63) 



9m 1 



9^2 
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and 



1 da d^Q 

V{ui,u^) := [(u3 + Z/)--(Mi + Z/)](t + — ^)-2(mi-U3)(u3 + Z^) ^ ^ =0. 

4 0U2 OU2OU3 

(64) 



Solving for t from (63) and substituting it into (64), we use (48) to simplify 
the result and get 

1 d'^Q 
W{ui,u-i) := [(% + z/) - -(Ml + z/)](mi + z/) . ^ 

4 OU1OU3 

+ [(Ml + z/) + 2(m3 + ^^)] (m3 + 1^) = • (65) 

OU2OU3 



Obviously, equations (63) and (64) are equivalent to equations (63) and (65) 



We now solve equation (65) for as a function of ui in the neighborhood of 



Ui = U3 = u. We use formula (50) and the symmetry of q to write 



d'^q{ui,U3,U3) 
duidus 

d'^q{ui,U3,U3) 
duodu-:!. 



1 



64v^ 



71 J-1 



1 „,l+/i ^ 



(66) 
(67) 



For Ml = M, it follows from (55), (66) and (67) that equation (65) has only the 



solution M3 = M in the neighborhood of Mi = M3 = u. 

For Ml which is a little bigger than m, we will show that there is a unique M3 



such that equation (65) holds. By (55), (66) and (67), we have 



Hence, 



W{u, M3) > for some M3 < m . 
Vr(Mi, M3) > for Ml a bit larger than u . 



For each of such Mi's, we deduct from (55) and (65) again that 



W{ui,u) < . 

By the mean value theorem, we show that, for each mi that is slightly larger 



than it, there exists a M3 < m such that (65) holds. It is easy to check the 
uniqueness of M3. 



Therefore, (65) determines M3 as a function of Mi, M3 = ^(mi), for small non- 
positive Ml with A{u) = u. The smoothness of the function /(m) and Lemma 
14 imply that A{ui) is a smooth decreasing function of Mi. 
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Next, substituting = A{ui) into (63), it is not hard to show that (63) 



determines Ui as a function of t. We have therefore proved the short time 
version of the following lemma. 



Lemma 15 Under the conditions of Theorem 13, equations (65) and (64) 
have a unique solution (u{(t), u^if), u^{t)) with M2 (t) = u^if) for all t > 0. 
The solution has the property that Ui {t) > U2 {t) = {t) for t > and that 
uTiO) =7x2(0) = ^3(0) =M. 



Proof. We will now extend the solution {ui (t), U2 (t), % (t)) of equations (63) 



and (64) for all t > 0. Before doing this, we need a lemma. 



Lemma 16 Under conditions of Theorem 13, the following hold: 



d'^q 



duidu2 duiduz 



d^q dq 

<0, 7ri<0, t+-^<0 
ouf OUi 



and 



t + 



dq 



dq 



at the solution (mi,^,^) of (65) where Ui > u^. 



Proof. By ( 48 ) and Lemma 14 



d'^q 



d'^q 



d dq dq 



duidus du2du3 du-s dui du2 



2(mi - U2] 



d^q 
duidu2du3 



< , (68) 



which when combined with ( 65 ) gives 



d'^q 



duiduz 



< 



as long as 4(^3 + z/) — {ui + u) > 0. This inequality holds even when 4(m3 + 
z/) — (ui + u) < 0. To see this, suppose the inequality fails at some point, at 



which d^ q must vanish because of (65). This would violate (68) 



The other inequalities of Lemma 16 are shown in the same way. 



□ 



We now calculate the partial derivatives of U and V at the solution {ui,Us, U3) 



of (63) and (64) 
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dU 
dui 
dU 

dU3 



dV 

du?, 



dq , , 9^0 , , S^q 

t+7T-^ + (Ml + Z^)7^+2(M3 + I/)^^ ^ <0, 



2(t + 



1 



dq 



Ui - u-s 
dq 



d'^q d^q 
) 2(ni + u) I + 8(m3 + i^) ^ I 
OU1OU3 OU2OU3 

{ui + iy){t+ p-) + 2(m3 + z/)(t + 



dui 



dm 



t+^ + [(8(n3 + z/)-3(ui + z/)]- 



d'^q 



du2du3 

d^q d^q 
-2K-.3)(% + -)(^^ + ^^) 

3(^1 + vf - 12(mi + z/)(m3 + v)+ 24(m3 + 

8(mi - M3)(m3 + l^) 

d^q d^q 



dq 
duo 



2{ui - M3)(m3 + 



duldu?. 



du2du\ ^ ^ ' 



where we have used Lemmas 14 and 16 to determine the signs of the deriva- 
tives. These show that the Jacobian 



d{Ui,U3) 



7^0 



on the solution {ui,Us,U3) of equations (63) and (64) where Ui > u^. 



Therefore, by the Imphcit Function Theorem, equations (63) and (64) can be 
solved for Ui(t), ^2 (t) = u^it) for all t > 0. Furthermore, it is easy to check 
that Ui{t) is an increasing function of t. □ 



5.2 The leading edge 



At the leading edge, Ui = U2, i.e., s = 1, it follows from (A. 7), (A. 8), (58), 

dq 



and (59) that equations (|60|) turn out to be 
dq 



t + ^iui,Ui,U3) 

dui 







dus 



. 



In the same way as we handle Lemma [15} we can solve the above equations 
for Ui and U3 as functions of t, leading to the lemma. 



Lemma 17 Under the conditions of Theorem 13, system (60) has a unique 
solution (ui{t), U2{t), u^lt)) with Ui{t) = U2{t) for all t > 0. The solution 
has the property that uf(t) = 'U^(t) > u^{t) for t > and that uf{0) = 
ut{0) = ut{0) = u. 
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5.3 Near the trailing edge 



By Lemma 15, {ui (t), U2 (t), (t)) satisfies equations (63) and (64). For each 



fixed t > 0, we need to solve equations (60) for ui and Us as functions of U2 



in tlie neigliborliood of U2 (t). This is carried out in 

Lemma 18 For each t > 0, equations (60) can be solved for Ui and U3 in 
terms of U2 in the neighborhood of (ui(t), ^2 (t), u^{t)) 

^ ' (69) 

M3 = N{U2) 



such that Ui (t) = M{u2 (t)) and % (t) = N{u2 (t)), Moreover, foru2 > U2 {t), 

N{U2) <U2< M{U2) . (70) 



Proof. Calculating the first partial derivatives of F and G/s of (61 ) and (62) 



at {ui (t), U2 (t), (t)), where U2 (t) = (t), and using (48), we find 



dF 

dui 
dF 



dF 



d{G/s^) 



djC/s' 
du?. 



d'^q _ 2{u2 + y) 
duf {ui + vY 
d\ , 2 

+ [: 



it 



dq ^ 2{u2 + v) d'^q 
du2 ui + u duidu2 
3 1 dq 



< 



duidu2 ^Ui + u Aui—Us 
6{u2 + v) d'^q 



du-2 



+ 



+ 



d'^q 



ui + u du2du3 







+ 



dq 2{u2 + y) d'^q 



+ 



duidus 4 ui — U3 du2 
TT j 19(ui + u) - 16(m3 + u 
2 j 32{ui - U3){ui + u) 
2{ui + vf + 9(Mi + z/)(m3 + z/) - 8(m3 + vf 



Ui + u du2du3 
dq ui - U2 
du2 2{ui + v) du2du3 



-0, 
d'^q 



TT 

2 



64(-ui - us){ui + z/)(m2 + ly) 
3 3(^1 + u) 



TT{t 



dq 
dui 



> 



16(^1-^3) 32{ui - us){u2 + ly)'"^" ' du2 
?>{u2 + u) d'^q [ui ~ u^){u2 + y) d^q 



2(mi + u) du2du^ 



Ui + V 



du2du% 



TT [ 3[(m + 



4(mi + z/)(m3 + + 8(^3 + uf] ^ dq 



32(mi - U'i){ui + h'){u2 + V 
{ui - u-i){u2 + y) d^q 



duo 



Ui + V 



du2dul 



> 
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where we have used (63) and (64) to simphfy the results, and Lemmas 14 and 
16 to determine the signs of the derivatives. 



These prove the non-vanishing of the Jacobian 

d{F,G/s^) 

d{ui,U3) 



at {ui (t), U2 if), U3 it)). Hence, equations (60) can be solved for 



Ui = M{U2) 
Us = N{U2) 



in a neighborhood of U2 it) such that Ui (t) = M{u2 it)) and ^3 (t) = N{u2 it)). 



Furthermore, N{u2) is a decreasing function of U2 and so (70) holds. □ 



5.4 The passage from the trailing edge to the leading edge 



We shall show that, for each fixed t > 0, solutions (69) of equations (60) can 



be further extended as long as N{u2) < U2 < M{u2). The Jacobian of system 



(56) and (57) with respect to (-Ui, ^3) has to be estimated along the extension. 



Lemma 19 Under the conditions of Theorem 13_, the following inequalities 
hold for each t > 0. 

9(Ait + «;i)^Q d{\2t + W2) ^ Q diM + w^^^ ^^^^ 



dui 



dui 



dm 



on the solution (ui, U2, u^) of (56) and (57) (or equivalently (60)) in the 



region < M3 + z/ < M2 + z/ < mi + i/. 



Proof. Using formulae (51) for wi, W2, and W3, we see that (56) and (57) are 
equivalent to 

[Al-(Ml+M2 + M3 + 2z/)](t+^) = [A2-(Ml+M2+M3 + 2z/)](t+|^), (72) 

[X2-2{ui + U2 + us + 2iy)]{t + p-) = [A3 - (ui +M2 + ^3 + 2z/)](t+ . (73) 

OU2 0U3 

By Lemma [161 



d^q{ui,U2,U3) d'^q{ui,U2,U3) ^ ^ 



duidu2 



duidu3 



(74) 
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at the trailing edge. 



We claim that inequalities (74) hold on the solution (■ui,M2,M3) of (56) and 
(57) as long as < + u < U2 + < Ui + u. 

We justify the claim by contradiction. Suppose otherwise, for instance at some 
point {ui,U2,U3) on the solution of (56) and (57), with 0<U3 + i'<U2 + i^< 

Ul + 1^, 



duidu2 



. 



In view of (48), this gives 



dq dq 



dui du2 



at (mi,M2,M3) 



which together with (19), (20), and (72) imply 

1 dq 



t + 



1 dq 
28^1 



t + 



2du2 



By ([20]), (12T|, ((73]), and (jTSj), we obtain 

dq 



t + 



which together with (48) gives 

d'^q 



d'^q 



duidu2 duidus 



at {Ui,U2,U3). 



On the other hand, by (48) and Lemma 14 



d'^q _ d^q ^ 2(^2 - u^) < 

duidu2 duidu^ duidu2du3 



(75) 



(76) 



at (mi,M2,M3). This contradicts (76) and the claim has been justified. 
By ( 48 ) , we have 

2(u - u ) dq 
duidus dui du3 ' 
Differentiating this with respect to ui yields 



d'^q 



d'^q 



dui duidu'- 



+ 2(mi - Us] 



d^q 
dujdus 



< 



(77) 



where we have used (74) and Lemma 14 in the last step. 
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It follows from ([74|) and (|48| that 



dq dq dq dq 



dui du2 ' dui du3 ' 



which when combined with (19), (20), (21), (72) and (73) gives 

^ 1 dq ^ ^ ^ 1 dq ^ ^ ^ 1 dq ^ ^ 
2 dui ' 2 du2 ' 2 (9^3 



(78) 



on the solution {ui, U2, M3) of (56) and (57) in the region < us + u < U2 + y < 

Ui + V. 



Therefore, by (51) 



' dui 



where in the last inequality we have used (19), (77), (78), and 



dui {duiiy 



> . 



This proves the first inequality of ( 71 ) . 



Next we shall prove the rest of Lemma 19, By (48), we have 

^ d^q dq dq 
du2dus du2 du3 ' 

Differentiating this with respect to U2 yields 

d'^Q o -92^ , . d^q 

^ = 3^-^ h 2{U2 - Us) . 

dU2 dU2dU3 dU2dU3 



(79) 



Using (48) to rewrite (73), we obtain 



da d'^a 
(A2 - A3)[t +^]+ 2[A3 - (Ml + U2 + U3 + 2z/)]— |-(m2 - %) = (80) 
dus du2du3 



which together with ( 79 ) gives 



0A2 - A3 . 5g d^q 
3 {t + — ) + 2[A2 - (Ml +U2 + U3 + 2z/)J^ 



U2 - U3 



dm 



4[A2 - (mi + U2 + U3 + 2z/)](m2 - Us) 



d^q 
du^dus 



> 



^1) 
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where we have used (21 ) and Lemma 14 in the last inequahty. 



It follows from (51) that 



d{\2t + W2) 0X2 



da 

it + ^) + [A2 - (Ml + M2 + M3 + 2Z/)] 



dul 



3 Ao — A.q 



> - 



dq 



2U2 - Us 
3 A2 — A3 

2 M2 - U3 

+ 3(A2-A3 



{t + 7^) + [A2 - (Ml +U2 + U3 + 2u)] 

Cm 2 
dq 

{t + ^) + [A2 - {ui +U2 + U3 + 2z/)] 



du3 
d'^q 



d'^q 
dul 
d^q 
dul 



du2dus 



>3(A2-A3 



d'^q 



du2dui 



> 



where we have used Lemma |4| (78) in the first inequality, and (81) in the 
second one. The last inequality is due to the fact that A2 — A3 and d^ q have 



the same sign because of (80) 



This proves the second inequality of (71). In the same way, we can prove the 
last one. □ 



We are ready to conclude the proof of Theorem 13 



Proof of Theorem 13 By Lemma 18, equations (60) can be solved for 



Ui = M{U2) 
Us = N{U2) 



in the neighborhood of {ui (t), U2 (t), M3 (t)). Furthermore, (70) holds if U2 



> U2 (t). We shall extend the solution in the positive U2 direction as far as 



possible. It follows from Corollary p] and Lemma 19 that, along the extension 



of (69) in the region ui + z/ > M2 + z/ > M3 + z/ > 0, the Jacobian matrix of 



(56) and (57) is diagonal and therefore is nonsingular. Furthermore, equations 



(56) and (57) determines (69) as two decreasing functions of U2- 



This immediately guarantees that (69) can be extended as far as necessary in 



the region -ui + z/ > M2 + z^ > M3 + z/ > 0. Since M{u2) is decreasing, (69 ) stops 
at some point 'U^(t) where, obviously, M{u2{t)) = U2{t). Therefore, we have 



shown that (56) and (57) determine Ui and ^3 as decreasing functions of U2 



over the interval [u2 (t),U2(t)]. 



28 



Let 



ufit) = M{ut{t)) 
utit) = Niutit)). 



Clearly, {ui (t) , U2 (t) , (t)) solves system (60) at the leading edge ui = U2 



Hence, these Mi (t), (^) ^-^d u^{t) are exactly the ones appearing in Lemma 



Substituting (69) into (44), we obtain 



X = A2(M(m2), M2, N{u2))t + W2{M{U2),U2, N {U2)) 

which by Corollary |2] and Lemma [19] clearly determines x as an increasing 
function of U2 over interval [u2{t),U2{t)]. It follows that, for each fixed t > 0, 
U2 is a function of x over the interval [x^{t)^ and that therefore so are 

Ml and M3, where 



x^ it) = A2 (nf it) , ut it) ,ut{t))t + W2 (nf (t) , (t) , (t) ) 



Thus, (44) can be solved for 



(82) 



Ui=Ui{x,t), U2 = U2{x,t), U3=U3{x,t) 

within a wedge 



X'(t) < X < x+(t) 
x-(0) =a;+(0) = 



for t > 



^3) 



where we have used (\82h, Lemma 15, and Lemma 17 in the last equations. 



Boundary conditions rt9| and ( 10 ) can be checked easily. The proof of Theorem 



13] would be completed if we can verify that the wedge is indeed a cusp. First 
we need a lemma. 



Lemma 20 At {ui (t), U2 (t), (t)), 

d{X2t + W2) _ diXst + W3) 
du2 du^ 

while at (uf(t), u^(t)), 

d{\it + wi) _ d{\2t + W2) 
du\ du2 



0. 



^5) 



Proof. 
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Using expansions (A. 5) and (A. 6), we obtain from (25) that 



dm 



at (mi (t), U2 it), % it)). 
By (fsTl), we have 



d{\2t + w2) ^ Q[{m + y) - \{:ui + p)] ^^ ^ dq_ 

{ui + u) dU2 
4(mi - M3)(m3 + v) d'^q 



Ui + u 







where we have used (22) in the first equahty, and (64) and (79) in the last 



equahty. The proof also applies to the other equation of (84) 



To prove (85), we proceed as follows. By (18), (A.7), and (A.8), we find 

8X2 



dui 



0{K{s)) as s ^ 1 



(86) 



By Corollary [2 



of (56) and (57) 



and formula (51), we calculate the derivative on the solution 



d{X2t + W2) 
dui 

8X2 , dq . 1 ,^ , ^ d'^q dq 



dq 



which when combined with (23) and (86) gives 



lim[t + 



dq 



K{s) = 0. 



On the other hand, as in (86), we have 



dX2 

duo 



0{K{s)) ass^l . 



Therefore, by (51) we see that at {ui, U2, u^) 



^7) 



d(X2t + W2) dX2 , 1 9o , / 



d'^q 
dui 
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where, by (87) and (88), the first term vanishes, while the second term vanishes 
because of (23). This proves the second equality of (85). In the same way, we 
can check the other equality of (85). □ 



Now we continue to finish the proof of Theorem 13 Differentiating (82) with 



respect to t, by Corollary [2] and Lemma 20 we obtain 



dx^t) d{X2t + W2) dui 



dt 



dU2 

A2(m^,m^, 



dt 



Un 



which when combined with (22), (23), Lemma 15, and Lemma 17 gives 



dx^{t) 
dt 



3m + 2z/ at t = . 



Therefore, wedge (83) is a cusp. This completes the proof of Theorem 13 



We immediately conclude from Theorem[T]and Theorem 13 the following result 
on the initial value problem of the Whitham equations. 



Theorem 21 For a decreasing initial function Uq{x) whose inverse function 
f{u) satisfying the conditions of Theorem 13, the Whitham equations have 
a solution {ui{x,t),U2{x,t),U3{x,t)) within a cusp for all positive time. The 
Burgers solution of ^ exists outside the cusp. The Whitham solution matches 
the Burgers solution on the boundary of the cusp in the fashion of and (10). 



We close this paper with two observations. First, it is obvious from the proof 



of Theorem 21 that one should expect local (in time) results if local conditions 

< for all 6 < M < a in 

< in the neighborhood 



are assumed. Namely, if the global condition /" 



Theorem 21 is replaced by a local condition /" 



of the breaking point u, the results of Theorem 21 are only true for a short 
time after the breaking time. 



Second, a hump-like initial function can be decomposed into a decreasing and 
an increasing parts. It is known that the decreasing part causes the Burgers 
solution of ^ to develop finite time singularities while the increasing part 
does not. These two pieces of data would not interact with each other for a 
short time after the breaking of the Burgers solution. As a consequence, a 
short time result also holds for a hump-like initial function. 
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A Complete Elliptic Integrals 



In this Appendix, we list some of tlie well-known properties of the complete 
elliptic integrals of the first, second and third kind. 



These are the derivative formulae 



dK{s) 


E{s) - (1 - s)K{s) 


ds 


2s(l - s) 


dE{s) 


E{s)-K{s) 


ds 


2s 


dli{p, s) 


pEis) + is-p)K{s) + ip'-s)U{p,s 


dp 


2p{l-p)ip-s) 


dn(p,s) 


E{s)-{l-s)U{p,s) 


ds 


2{1 -s){s-p) 



(Al) 
(A2) 
(A.3) 
(A.4) 



K{s) and E{s) have the expansions 



i?(.) = "[l-^- ^( ^•^•••('"-^) )V----] (A.6) 

^ ^ 2 ^ 4 64 2n-V 2 ■ 4 ■ ■ ■ 2n ^ ^ ^ ^ 

for \s\ < 1. They also have the asymptotics 



K{s) ^hog^^, (A 7) 

E{s)^l + \{l-s)[log^^-l] (A.8) 

when s is close to 1. They satisfy the inequalities [H] 

1 K(s) 1-1 , , , , 

< < /or < s < 1 . (A.9) 



1-f E{s) 1-s 
The complete elliptic integral of the third kind has the following behavior 



n(p,s) = | whenp = 0, s = 0, (A.IO) 

— ^^7^ ~ when s is close to 1 . (A. 11) 

K[s) 1 — p 
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